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Preface 


One  expects  to  *ake  discoveries  as  a  result  of  research.  In  this  project  I 
discovered  that  a  thesis  is  not  something  which  sprouts  like  an  acorn  and  grows 
slowly  to  maturity.  This  one  grew  more  like  a  bacterial  culture:  ideas  flourished, 
only  to  be  proved  wrong  or  irrelevant  and  be  wiped  out  with  electronic  alacrity. 
Countless  glorious  inspirations  have  met  their  doom  since  I  naively  wrote  my 
prospectus;  in  fact,  I  imagine  no  more  than  ten  percent  of  even  my  "final"  thesis 
draft  is  still  included.  It  was  fun,  but  I  hope  to  imbue  my  next  project  with  a 
little  less  of  that  sort  of  excitement. 

The  excitement  could  have  been  a  killer  if  I  hadn't  had  a  lot  of  help.  My 
thesis  advisor,  Lieutenant-Colonel  Jim  Lupo,  saved  me  a  lot  of  frustration.  When 
I  bogged  down,  he  always  had  an  answer.  The  other  members  of  my  committee. 
Captain  Haaland  and  Major  Beller,  performed  the  invaluable  if  painful  service  of 
bringing  me  face  to  face  with  my  many  errors  and  lapses.  Captain  Steve  Thompson 
of  the  Air  Force  Astronautics  Laboratory  was  an  indispensable  source  of  information. 
I  am  also  grateful  to  Captain  Frank  Turner  for  sharing  his  knowledge  of  things 
academic  and  otherwise. 
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Abstract 


'A  computer  model  of  a  cryogenic  system  for  storing  solid  antimatter  is  used 
to  explore  the  radiative  cooling-power  requirements  for  long-term  antimatter 
storage.  If  vacuum-chamber  pressures  as  low  as  torr  can  be  reached,  and 

the  rest  of  the  large  set  of  assumptions  is  valid,  milligram  quantities  of  solid 
antimatter  could  be  stored  indefinitely  at  1.5  K  using  cooling  powers  of  less  than 
a  microwatt.  Many  of  the  assumptions  made  are  problematic  and  need  verification, 
as  they  could  potentially  change  the  results  greatly.  The  system  modeled  is  a 
sphere  of  solid  anti-parahydrogen  at  1.5  K  or  below  levitated  in  a  spherical 
cryogenic  vacuum  chamber.  The  free  matter  gas  in  the  chamber  is  assumed  to  be 
molecular  hydrogen,  and  sublimation  of  both  matter  and  antimatter  is  assumed  to 
be  negligible.  The  antihydrogen  is  assumed  to  be  in  thermal  equilibrium,  although 
annihilation-energy  deposition  is  localized  and  hydrogen's  thermal-impulse  response 
time  is  comparable  to  the  interval  between  annihilations.  A  parametric  analysis 
is  performed,  with  system  cooling  power  evaluated  over  a  wide  range  of  pressures 
and  system  sizes,  as  temperature  and  emissivity  are  also  varied.  The  cooling-power 


requirements  of  storing  solid  antimatter  for  extended  times  may  not  be  an  obstacle, 
if  the  proper  conditions'  o§tairf^  However,  whether  these  conditions  are  indeed 
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I  Introduction 


This  project  is  a  study  of  the  cooling-power  requirements  of  a  storage  system 
for  solid  antihydrogen.  The  approach  is  analytical,  beginning  with  the  derivation 
of  the  differential  equations  governing  the  system's  cooling  power.  The  equations 
are  simplified  and  solved  under  quasi-steady-state  conditions,  and  the  results 
incorporated  into  a  computer  program  that  evaluates  cooling  power  requirements 
as  various  system  parameters  are  varied. 

1  Cryogenic  Storage  of  Antimatter 

If  a  number  of  technical  obstacles  can  be  overcome,  the  matter-antimatter 
annihilation  reaction  may  be  useful  as  an  energy  source  for  applications  requiring 
compact  fuel  supplies.  For  example,  no  other  combination  of  potential  rocket 
propellants  releases  more  energy  per  unit  mass  of  fuel.  One  of  the  many  problems 
which  must  be  solved  if  antimatter  is  to  become  a  practical  fuel  is  how  to  store 
antimatter  so  that  it  does  not  react  with  matter.  A  very  high  vacuum  and  very 
low  temperatures  are  required  to  maintain  control  of  the  antimatter  in  storage. 

The  critical  component  of  antimatter  fuel  is  the  antiproton.  Unlike  the  more 
easily-obtained  positron,  the  antiproton  produces  energetic  charged  particles  when 
it  annihilates.  The  momentum  of  these  particles  can  interact  with  a  magnetic  field 
to  accelerate  a  rocket.1  Stored  antiprotons  must  be  prevented  from  touching 
normal  matter  while  in  storage,  or  they  will  annihilate  prematurely.  The  typical 
proposed  way  to  maintain  control  of  an  aggregate  of  antiprotons  is  to  combine 
them  with  positrons  to  produce  antihydrogen  atoms  and  create  a  pellet  of  solid 
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molecular  antihydrogen  in  a  vacuum  chamber.  It  may  not  be  possible  to  solidify 
a  gas  without  a  nucleation  site  as  would  be  necessary  with  antimatter;  research 
into  this  problem  is  under  way. 2  if  a  solid  antihydrogen  pellet  can  be  created, 
the  ice  ball  can  be  placed  in  a  cryogenic  vacuum  chamber  and  levitated  elec¬ 
trostatically  or  magnetically  to  keep  it  from  touching  the  chamber  walls. 
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Figure  1-1.  Vapor  Pressure  of  Solid  Parahydrogen 


The  ice  ball  must  be  maintained  at  a  very  low  temperature  in  order  to  keep  the 
antihydrogen  from  subliming.  At  temperatures  below  1  K  the  vapor  pressure  of 
hydrogen  is  negligible,  so  there  is  no  sublimation  and  the  antimatter  is  effectively 
contained.  Figure  1-1  shows  an  estimate  of  the  behavior  of  the  vapor  pressure 
between  1  and  3  K.  derived  from  a  cubic  least-squares  fit  to  published  calculations 


2 


of  the  pressures  at  1,  2,  3,  and  4  K.3  As  the  ice  ball's  temperature  exceeds  1  K, 
the  vapor  pressure  increases  rapidly,  prompting  increased  sublimation  rates  and 
eventually  allowing  significant  amounts  of  antimatter  to  reach  the  wall  and 
annihilate.  Some  of  the  annihilation  radiation  is  absorbed  by  the  ice  ball,  raising 
its  temperat  :e  and  causing  further  increases  in  sublimation  until  the  system  fails 
catastrophically. 

The  very  low  temperatures  and  pressures  required  for  this  method  of  storing 
antimatter  may  be  attainable  with  current  technology.  Pomeranchuk  cooling  systems 
using  3He  are  capable  of  sustaining  1  microwatt  of  cooling  power  at  25  mK.4  At 
such  low  temperatures,  the  nearly-perfect  vacuum  required  becomes  theoretically 
possible.  Sublimation  and  outgassing  of  solids  are  nil,  and  any  gas  particles  free 
in  the  chamber  will  adhere  to  the  walls  on  contact.  Experiment  has  shown  unpumped 
vacuum  chamber  pressures  falling  rapidly  below  10“ 14  torr  (the  minimum  measurable) 
as  temperature  decreases  below  30  K.5 

2  Scope  of  Study 

The  storage  system  modeled  in  this  study  is  a  spherical  cryogenic  vacuum 
chamber  containing  a  spherical  antihydrogen  ice  particle.  Several  system  parameters 
are  varied  and  their  effects  on  cooling  power  are  numerically  computed.  The  origin 
of  the  antihydrogen  sphere  and  its  removal  from  the  chamber  are  not  addressed. 
Temperatures  of  the  antihydrogen  ice  are  allowed  to  vary  from  millikelvins  to  1.5 
K.  If  the  temperature  of  the  antihydrogen  rises  above  1.5  K,  the  system  is 
considered  to  have  failed;  no  analysis  is  done  of  the  behavior  of  the  system  when 
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sublimation  becomes  significant.  The  following  assumptions  are  made: 

1.  The  ice  ball  and  the  chamber  are  each  characterized  by  a  single  temperature. 
The  form  of  the  radiative  heat  transfer  equation  used  is  invalid  if  temperatures 
are  not  uniform.  This  assumption  is  correct  if  the  smallest  time  intervals 
considered  are  much  larger  than  the  time  needed  to  conduct  heat  through 
the  chamber  wall  and  the  ice  ball.  The  thermal  equilibration  time  of  a  0.6-cm 
hydrogen  sphere  is  on  the  order  of  1  ms. 

2.  The  ice  ball  is  pure  anti-parahydrogen.  This  is  the  equilibrium  condition 
at  the  temperatures  considered  in  this  study,  but  the  nearness  to  this 
condition  in  a  real  system  will  depend  on  the  method  used  to  create  the 
antihydrogen. 

3.  The  chamber  wall  can  be  modeled  as  pure  iron  for  purposes  of 
ionizing-radiation  absorption  analysis.  A  typical  material  for  vacuum  chambers 
is  type  304  stainless  steel,  which  contains  about  35  parts  per  million  of 
nonferrous  elements.4 5 6 

4.  The  chamber  wall  is  coated  on  the  Inside  with  at  least  a  molecular 
monolayer  of  hydrogen,  so  that  its  radiative  graybody  emissivity  is  the  same 
as  that  of  hydrogen. 

5.  All  of  the  matter  particles  free  in  the  vacuum  chamber  are  hydrogen 

molecules.  Thi  would  be  approximately  true  if  the  vacuum  system  had  been 
pumped  out  sufficiently. 
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6.  Through  some  leakage  or  diffusion  process,  the  chamber  remains  continuously 
filled  with  hydrogen  gas  at  constant  pressure  and  at  the  temperature  of  the 
chamber  wall. 

7.  All  collisions  between  hydrogen  and  antihydrogen  molecules  result  in 
immediate  annihilation  of  both  proton-antiproton  pairs  and  both  electron- 
positron  pairs.  There  is  no  scattering  at  the  low  interaction  velocities 
characteristic  of  a  1-K  environment. 

8.  Both  the  ice  ball  and  the  chamber  wall  are  "well-behaved"  radiative 
graybodies  with  diffuse  reflection. 
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II  Computation  of  Cooling  Power  Requirements 

To  hold  the  vacuum-chamber  wall  at  some  fixed  temperature  T any  heat 
energy  absorbed  by  the  wall  must  be  removed  at  the  same  rate  it  is  deposited. 
The  chamber  wall  can  absorb  energy  from  three  sources:  heat  radiation  from  the 
ice  ball  in  the  center  of  the  chamber,  annihilation  radiation  from  matter-antimatter 
reactions  on  the  ice  ball  and  on  the  wall,  and  adsorption  and  reverse  sublimation 
of  vacuum  particles  onto  the  wall.  Besides  the  energy  removed  from  the  wall  by 
the  cryogenic  system,  the  wall  can  also  lose  energy  by  sublimation  and  desorption 
into  the  chamber. 

The  extreme  low  temperature  of  the  system  under  consideration  allows  con¬ 
siderable  simplification.  At  such  temperatures,  there  is  no  sublimation  or  desorption. 
Since  adsorption  of  particles  still  occurs,  the  vacuum  pressure  can  become  arbitrarily 
low  as  gas  particles  adhere  to  the  chamber  wall  and  are  not  replaced.  Since  the 
antihydrogen  does  not  sublime,  no  antimatter  is  free  to  annihilate  on  the  chamber 
wall.  Therefore  the  central  sphere  is  the  only  source  of  annihilation  radiation, 
and  this  source  only  lasts  as  long  as  it  takes  for  all  residual  particles  to  hit 
either  the  chamber  wall  or  the  sphere.  For  this  study,  the  system  is  assumed  to 
contain  some  hydrogen  gas  at  a  constant  but  very  low  pressure. 

1  System  Cooling  Power 

The  required  system  cooling  power  can  be  derived  from  an  energy  balance 
relation  for  the  chamber  wall.  The  chamber  wall  temperature  T2  is  assumed  to  be 
constant.  Therefore, 
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dT  2 
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[/3a2-/>s2  +  Q-/5r]=0 


(1) 


where  is  the  mass  of  the  chamber,  cV2  is  the  specific  heat  of  the  chamber 
material,  Pa2  is  the  power  deposited  in  the  wall  by  annihilations  on  the  inner 
sphere,  PS2  is  power  removed  by  sublimation  and  desorption  from  the  wall,  q  is  the 
rate  of  heat  energy  transfer  from  the  inner  to  the  outer  sphere,  and  Pf  is  the 
total  cooling  power  requirement.  Because  of  the  low  temperature  of  the  system, 
there  is  no  sublimation  or  desorption,  however.  Solving  the  above  equation  for 
Pj  (cancelling  out  ®nd  Cv2>  an<i  neglecting  the  sublimation  term. 


PT  =  P 


a  2 


+  q 


(2) 


Assuming  both  surfaces  to  be  gray  and  reflection  to  be  diffuse,  the  rate  of 
heat  energy  transfer  from  the  inner  to  the  outer  sphere  is7 


4nr^o (T,  -  T\) 


where  a  is  the  Stefan-Boltzmann  constant,  Tj,  Ap  r and  e ,  are  the  temperature, 
surface  area,  radius,  and  emissivity  of  the  ice  ball,  and  the  corresponding  quantities 
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subscripted  with  a  2  refer  to  the  chamber  wall.  Strictly  speaking,  the  substances' 
emissivities  depend  on  temperature,  but  over  the  narrow  temperature  range  under 
consideration  here  it  is  justifiable  to  assume  they  are  constant. 

The  power  deposited  in  the  wall  by  annihilations  in  the  central  sphere  is 
given  by 

2  =  (£a,*£„2)/?„,  <4> 

where  Raj  is  the  annihilation  rate  in  the  central  sphere,  Egj  is  the  energy 
deposited  in  the  central  sphere  per  annihilation  in  the  sphere,  and  £a2  is  the 
energy  deposited  in  the  wall  per  annihilation  in  the  sphere. 

The  annihilation  rate  in  the  inner  sphere  is  simply  the  collision  rate  for 
the  vacuum  gas  with  the  surface  of  the  inner  sphere,  derivable  from  the  kinetic 
theory  of  gases:8 


Here  nv  is  the  vacuum  number  density,  vm  is  the  mean  gas  particle  velocity,  Pv  is 
the  vacuum  chamber  pressure,  mv  is  the  mass  of  a  gas  molecule,  Tv  is  the  gas 
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temperature,  and  k  is  Boltzmann's  constant.  At  a  fixed  chamber  pressure,  there  is 
thus  an  inverse  dependence  of  the  collision  rate  on  gas  temperature,  devolving 
from  the  ideal-gas  law  relating  temperature,  pressure,  and  number  density. 

2  Annihilation  Energy  Deposition 

The  mutual  annihilation  of  a  hydrogen  atom  and  an  antihydrogen  atom  liberates 
the  rest  mass  energy  of  two  protons  and  two  electrons,  or  1877  MeV.  Since  the 
free  particles  in  the  vacuum  chamber  are  assumed  to  be  diatomic  hydrogen  molecules, 
each  collision  yields  3755  MeV.  The  energy  from  a  proton-antiproton  annihilation 
appears  in  the  form  of  charged  and  neutral  pions.  Forward®  reports  that  each 
proton-antiproton  annihilation  releases  its  energy  on  the  average  in  the  form  of 
3.0  charged  pions  and  1.5  neutral  pions,  (as  well  as  negligible  amounts  of  other 
products:  0.08  kaons  and  0.02  prompt  gamma  rays  per  annihilation).  The  neutral 
pions  decay  immediately  to  high-energy  gamma  photon  pairs.  The  energy  from  an 
electron-positron  annihilation  appears  as  two  (or  rarely  in  solids,  three)  gamma 
photons.  The  energy  deposition  modes  to  be  considered  are  thus  charged  pions, 
pion  decay  gammas,  and  annihilation  gammas. 

lie  Antihydrogen  Sphere 

Charged  Pions 

Figure  II- 1  shows  the  energy  spectrum  of  the  charged  pions  from  a  low-energy 
proton-antiproton  collision.  Gaines1®  divides  the  pion  spectrum  into  approximately 
70  energy  groups  and  finds  the  average  stopping  power  of  solid  hydrogen  for  a 
single  charged  pion  from  this  energy  spectrum: 
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Probability 


Energy  (MeV) 


- Charged  Pions  —  Pion-Decay  Gammas 

Mean  Energy  2S0  MeV  Mean  Energy  200  MeV 

Figure  II- 1.  Cbarged-Pion  and  Gamma  Energy  Spectra 


0.437  MeV/cm 


-  7.00  •  10~12  J/m. 


(6) 


This  formulation  assumes  the  plon's  energy,  and  hence  the  stopping  power,  remains 
nearly  constant  as  the  particle  traverses  the  hydrogen.  This  can  be  seen  to  be 
valid  for  the  pion  energies  (on  the  order  of  hundreds  of  MeV)  and  hydrogen  path 
lengths  (on  the  order  of  millimeters)  under  consideration  here. 
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Given  the  average  energy  deposition  per  pion,  the  i»ean  energy  deposited  per 
molecular  annihilation  can  be  calculated.  Each  molecular  annihilation  produces  an 
average  of  six  charged  pions,  three  from  each  proton-antiproton  pair.  Since  the 
reaction  products  are  emitted  isotropically,  only  half  the  total  number,  or  three 
per  molecular  annihilation,  pass  through  the  hydrogen  sphere.  Thus  the  energy 
deposited  in  the  hydrogen  by  a  single  molecular  collision  is 

fR1  -3.00- (7.00-  10'l2J/m)<x>  (7) 

-(2.10-  10~M  J/m)-  - 

rt 

where  d  is  the  diameter  of  the  sphere  and  <x>  -  d/n  is  a  mean  chord  length 
through  the  sphere.  This  analysis  differs  from  Gaines'  in  using  <x>  instead  of 
the  full  diameter,  and  in  assuming  molecular  rather  than  single-proton  collisions. 

Gamma  Rays  from  Neutral-Plon  Decay 

The  neutral  pions  decay  to  pairs  of  gamma  rays  after  a  mean  lifetime  of  just 
9.0  x  10~17  seconds,  before  they  can  deposit  any  significant  amount  of  energy  in 
the  system.  The  gamma  rays  resulting  from  the  pion  decay,  however,  do  deposit 
energy.  The  energy  spectrum  of  the  pion-decay  gammas  is  also  given  in  Figure 
II— 1.  Gaines  uses  the  decay-gamma  spectrum  in  Figure  II-l  to  compute  an  average 
linear  gamma  attenuation  factor  of  0.198  MeV/cm.  Thus  for  a  molecular  annihilation, 
three  neutral  pions  will  be  produced  which  will  decay  to  six  gammas,  three  of 
which  will  penetrate  the  hydrogen  and  deposit 
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£YRl  =3.00-  (0.1 98MeV/cm)-  <x> 


(8) 


=  (9.52-  10~l2J/m)  • 

n 


Figure  11-2.  Gamma  Attenuation  Coefficients  in  Iron  and  Hydrogen 

The  gamma-ray  portion  of  Gaines'  analysis  is  for  single-proton  collisions  with  the 
gammas  traversing  the  entire  diameter  of  the  sphere  as  before,  and  it  erroneously 
stipulates  that  each  neutral  pion  decays  into  one  gamma,  not  two.  Gaines'  linear 
model  of  gamma  energy  deposition  is  valid  because  the  size  of  the  hydrogen  ball 
is  very  saall  cospared  to  the  interaction  lengths  of  these  high-energy  gammas, 
as  shown  in  Figure  II-2.  For  example,  the  interaction  length  of  a  comparatively 
less-penetrating  10-MeV  gamma  is 


(9) 


=  (0.21  cm2/ g-  0.0893 g/cm3)'1 
1 

0.0  18cm 
»  53cm 

where  p- 0.0893  g/cm3  Js  the  density  of  solid  hydrogen  at  zero  pressure. 

Gamma  Rays  from  Electron-Positron  Annihilation 

In  solids,  the  annihilation  of  electrons  and  positrons  almost  always  produces 
a  single  pair  of  511-keV  gamma  photons.  The  energy  deposited  by  such  photons 
in  a  mean  chord  <x>  of  a  solid-hydrogen  sphere  is 


r  f  rn  ,  .  .  .  M*  /p/1  -<tl,  ',p)(p)<*>.. 

Ey.\  =(51  1  keV)-j^—  (1  -e  ) 


=  (8.19-  10~'4j)^YU-e 


-(0.1  2cm  */j  )  (0.08939 /cm  3)  (d/n) 


(4.09*  10"I4J)(1 -e'(0'0107cm’  )(d/rt)) 


do) 


where  the  factor  n„/p,  accounts  for  non-local  energy  deposition,  and  the  total 
and  energy-deposition  cross  sections  are  those  of  molecular  hydrogen.11  For 
reasonably  large  values  of  d,  the  energy  deposition  from  electron-positron 
annihilation  is  sufficiently  small  compared  to  the  energy  deposition  from  charged 
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pions  to  be  neglected.  For  example,  a  small  Ice  particle  of  2.0  ng  mass,  with  a 
diameter  of  0.0036  cm,  will  absorb  1.5  x  10“20  J  from  annihilation  photons  and  9.7 
x  10-14  J  from  charged  pions. 

Total  Annihilation  Energy  Deposited  in  Antihydrogen 

The  contributions  from  charged  and  neutral  pions  can  now  be  summed  (neglecting 
the  annihilation  gammas),  giving 

£,=£,  +  £  ,  (11) 
^  a  I  *-  n  l  L  ynl 

J 

=  (2. 1 0  •  10' 11  J/m  +  9.5-  1 0  ~ 12  J/m  )- 

ri 

=  (3.05-  10'"  J/m)-. 

ji 

as  the  expression  for  the  total  annihilation  energy  deposited  in  the  antihydrogen 
sphere  per  molecular  collision. 

Energy  Deposition  in  the  Chamber  Wall 

Gaines  does  not  extend  his  analysis  of  annihilation-energy  deposition  in 
hydrogen  to  the  deposition  of  this  energy  in  a  steel  wall.  The  same  techniques 
are  applicable,  but  since  the  wall  thickness  is  not  as  much  less  than  the  interaction 
lengths  for  pions  and  gammas,  more  accurate  alternative  methods  are  used  for  this 
study.  Annihilation  energy  strikes  the  wall  in  the  form  of  charged  pions  and 
pion-decay  gammas  from  annihilations  on  the  antihydrogen  sphere.  As  in  the 
antihydrogen,  gammas  from  electron-positron  annihilation  are  not  important  sources 
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of  energy.  The  total  annihilation  energy  deposited  in  the  wall  is 


E  ,  ■  E  o  +  E  o  (12) 

where  the  energy  deposited  by  charged  pions  is  determined  from  tables  of  pion 
ranges  in  iron,  given  an  assumed  thickness  x  of  the  chamber  wall. 


Charged  Pions 

In  a  residual -range  calculation,  the  energy  deposited  is  determined  by  finding 
the  range  of  the  pion  in  the  material  at  the  incident  energy  and  subtracting  the 
thickness  of  the  material  to  find  the  pion's  remaining  range  in  the  material.  The 
pion  energy  associated  with  this  range  is  found  from  the  range  table  and  stipulated 
to  be  the  pirn's  energy  upon  exiting  the  wall.  The  energy  deposited  is  just  the 
difference  between  the  entering  and  exiting  pion  energies.  This  process  is  repeated 
over  many  pion -energy  bands,  and  the  fractional  results  are  summed  to  reach  the 
total  pion  energy  deposited. 


Gamma  Rays  from  Neutral-Pion  Decay 

The  energy  deposited  by  pion-decay  gammas  is  also  calculated  by  summing 
the  contribution  of  many  energy  bands,  as 


i  M*  t 


(13) 


where  f/  is  the  fraction  of  the  spectrum  modeled  as  having  energy  £2\  Here  x 
represents  the  thickness  of  iron  traversed  by  the  radiation,  which  is  approximately 
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the  thickness  of  the  chamber  wall  since  the  radiation  originates  near  the  center 
of  a  spherical  chamber.  The  factor  p,„/pt  is  to  account  for  a  fraction  of  the 
energy  released  from  a  photon-matter  interaction  not  being  deposited  locally. 

This  device,  the  energy-transfer  ratio,  is  not  useful  at  the  high  photon 
energies  (greater  than  10  MeV)  under  consideration  here.  At  such  energies,  pair 
production  is  the  most  likely  interaction  mode.  Although  only  a  fraction  of  the 
incident  photons  will  interact  with  the  target  matter  along  a  given  linear  path 
length,  those  causing  pair  productions  may  transfer  nearly  all  of  their  energy  to 
the  resulting  energetic  electron-positron  pair.  The  energetic  electron  thus  produced 
may  deposit  all  its  energy  in  the  material  or  may  escape,  carrying  off  some  energy. 
The  positron  also  may  escape,  and  if  it  is  stopped,  it  will  be  annihilated  and 
produce  511-keV  photons,  which  may  also  escape.  Both  the  electron  and  the  positron 
emit  Bremsstrahlung  radiation  which  may  also  escape.  The  fraction  of  interaction 
energy  which  remains  in  the  system  is  thus  very  sensitive  to  the  materials  and 
geometry  of  the  wall.  As  a  worst-case  simplification,  the  ratio  p,„/p.,  is  assumed 
to  be  unity  in  this  study. 

3  Temperature  of  the  Central  Sphere 

The  inner  sphere's  temperatrre  is  governed  by  radiative  exchange  with  the 
chamber  wall,  and  by  the  annihilation -energy  deposition  rate.  If  the  ice  ball's 
mass  mj  is  assumed  to  be  constant,  the  rate  of  change  of  the  ice  ball's  temperature 
is 


16 


dr, 

d  t 


m  i  c, 


(Pa\~q) 


(14) 


where  Paj  represents  the  annihilation-energy  deposition  rate  on  the  inner  sphere. 
The  specific  heat  cvj  can  be  calculated  from  a  simplified  expression  of  the  Debye 
formula*2 


12n4S/'£V  1151 

C"'~  5 Mh  UJ 

since  the  temperatures  under  consideration  are  low  compared  to  hydrogen's  Debye 
temperature  0O-  120.3K.  The  factors  R  and  Mfj  in  the  Debye  equation  represent 
the  gas  constant  and  the  molecular  weight  of  hydrogen. 


4  Antihydrogen  Sublimation 

Equilibrium  Sublimation  Rate 

For  this  study,  a  temperature  range  is  chosen  within  which  there  is  effectively 
no  sublimation  of  solid  parahydrogen.  The  sublimation  rate  at  a  given  temperature 
can  be  estimated  in  two  ways,  which  together  provide  an  upper  and  lower  bound  on 
the  sublimation  rate. 


An  upper  bound  can  be  placed  on  the  sublimation  rate  if  one  assumes  that 
the  lattice  is  vibrating  at  its  maximum  frequency,  that  associated  with  the  Debye 
temperature.  In  this  case,  the  Arrhenius  relation  gives  the  resulting  sublimation 
rate  as*3 
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»  i 

kQAl 

'  rn 

kT  ,/k.T  . 

4jx  — 

e  *  1 

h  )\ 

,  v) 

4/c  n 

120. 3K  1 

h  (3.7-  1 0'lom  )2  J  1 


(  -  89.8  K  )/  7 


(16) 


where  v  is  the  vibration  frequency,  n*  is  the  number  of  molecules  on  the  surface, 
AE  ,is  the  enthalpy  of  sublimation,  h  is  Planck's  constant,  and  \s  is  the  intermolecular 
spacing.  At  a  given  temperature,  the  sublimation  rate  can  be  no  higher  than  that 
given  by  this  expression,  since  the  lattice  vibrations  can  be  no  faster  than  vm. 


The  sublimation  rate  can  be  bounded  below  by  reasoning  from  a  hypothetical 
case  where  the  hydrogen  or  antihydrogen  is  in  equilibrium  with  its  vapor.  Under 
such  conditions,  the  rate  of  sublimation  must  equal  the  rate  of  collision  of  vapor 
particles  with  the  surface.  The  sublimation  rate  will  therefore  be 


R 


4rtr | Pv 
V 2nmHkT , 


(17) 


where  Py  is  the  vapor  pressure  from  Figure  1-1  and  is  the  mass  of  a  hydrogen 
molecule.  The  solid  in  vacuum  should  not  sublime  any  more  slowly  than  it  does  in 
equilibrium. 
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Sublimation  Kate  (per  second) 


Temperature  (K) 

Equilibrium  Model  — Arrhenius  Model 

1  sublimation/year  —  1  sublimation/hour 

Figure  11-3.  Sublimation  Rates  in  a  l-mg  Solid- Hydrogen  Sphere 

Figure  II- 3  shows  the  two  estimates  of  the  sublimation  rate.  Based  on  these 
calculations,  it  seems  credible  to  consider  the  sublimation  rate  to  be  effectively 
zero  at  temperatures  as  high  as  1.5  K.  This,  however,  ignores  the  possibility  of 
localized  sublimation  caused  by  the  deposition  of  annihilation  energy. 

Local  Sublimation 

Local  energy  deposition  by  annihilation  radiation  may  be  sufficient  to  cause 
sublimation.  This  possibility  can  be  examined  using  the  stopping  power  and  linear 
attenuation  factor  computed  by  Gaines.  The  energy  required  to  sublime  a  single 
hydrogen  molecule  is14 
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(18) 


AEs  =  kTs 

=  (  1 .38  •  1  0  23  J/K )  (89.8K  ) 

=  1.24-  10'2' J 
=  7.74-  1 0"3eV . 

If  a  single  charged  pion  deposits  0.437  MeV  per  centimeter  of  hydrogen  traversed, 
and  the  intermolecular  spacing  is  3.7  angstroms,  then  each  molecule  along  the 
pion's  path  absorbs 


=  (4.37-  1 05eV /cm )  (3.7  •  10'8cm) 

-  1. 62-  1  0’2eV 
=  2.09-  A Es. 


Thus  if  all  of  the  absorbed  energy  went  into  heating  a  single  string  of  molecules, 
it  would  take  only  one  pion  track  to  cause  considerable  sublimation. 

In  fact,  however,  the  energy  is  not  deposited  along  a  one -dimensional  string. 
The  gamma  rays  from  neutral-pion  decay,  which  have  approximately  half  the  charged 
pions'  linear  energy  transfer  rate,  deposit  their  energy  throughout  the  solid.  A 
charged  pion  deposits  its  energy  in  a  cylindrical  annulus  of  dimensions  determined 
by  the  minimum  and  maximum  possible  impact  parameters  for  electromagnetic 
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the  minimum  impact  parameter  is  effectively  zero: 

ft  (21) 

bm'°-yrn.v 

he 

ymtc2\l  1  -  y'2 

1.973-  10~l3MeV  ■  m 
(2.79) (0.51  1  MeV  )>/  1  -  2.79  2 

=  1.48-  1 0"  l3m  . 


The  maxiaum  impact  parameter  is 
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ft  try  (22) 


ft  rr  ~ 

■"  1  "  V 

(6.582-  10'1  aV-  s)(2.79)(2.998  •  108rh/G)  n —  -~i 

- * : - — - — - V  1  -  2.79 

1  8.2e  V 

-  2.82  •  10'8m 


where  I  =  18.2  eV  is  the  mean  ionization  potential  of  molecular  hydrogen.16 
Considering  the  energy-deposition  region  as  a  cylinder  of  radius  bmax,  the  mean 
energy  deposited  per  molecule  is 


A£n=  (1.62-  10 


(3.7*  10~lom)2 

®  ’  )  ft  2 

4ji(2.82  •  10'8m) 


(23) 


=  2.21  •  1 0"7eV  . 

Averaged  over  the  entire  charged-pion  spectrum,  the  mean  energy  deposition  is 
4.71  x  10~7  eV  per  molecule,  still  well  below  the  energy  needed  for  sublimation. 
However,  since  energy  deposition  is  not  constant  throughout  the  cylinder,  but 
rather  varies  as  l/r^,  this  is  not  necessarily  a  conclusive  demonstration  that 
annihilation  radiation  will  not  stimulate  local  sublimation  of  the  antihydrogen. 


5  Heat  Transport  in  Solid  Hydrogen 

The  speed  with  which  the  added  heat  from  a  localized  energy  deposition 
travels  bears  on  the  validity  of  a  single-temperature  model  of  the  ice  ball.  In 
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order  for  the  model  to  be  valid,  the  time  required  for  the  heat  impulse  to  travel 


across  the  ball  must  be  small  compared  to  the  time  interval  between  energy 
depositions. 

The  ice  ball's  response  to  a  heat  impulse  is  given  by  the  Green's  function 
for  the  heat  equation:  1 7 

j\  f  a  —  -  n/2  -r2/4at  (24) 

G(r,0=(4na0  e 

where  a  is  the  thermal  diffusivity  of  the  material  and  n  is  the  number  of  dimensions. 
If  the  annihilation  impulse  is  deposited  along  a  narrow  track,  the  system  approximates 
a  two-dimensional  geometry,  so  that 


G(r,0  = 


4nA: ,  t 
P,c„l 


a/2 

-''2Pl  c„,/4k,t 

G 


4n(5M yrt/:  -r2p,(l2n**T3)/4(5M„e2)f 

(12n4/?r3)p  J  6 


(25) 


where  the  thermal  conductivity  kj  =  2.5  W/m/K  at  1.5  K.18  Figure  II— 4  is  a  graph 
of  the  Green's  function  at  T  =  1.5  K  and  r  =  0.6  and  1.2  cm,  the  radius  and  diameter 
of  a  10-mg  hydrogen  sphere.  The  figure  shows  the  impulse  response  peaking  at 
times  near  1  microsecond  after  the  impulse,  and  decaying  to  10X  of  its  peak  value 
after  less  than  100  microseconds.  For  a  10-mg  sphere  at  1.5  K,  a  collision  rate 
of  1  per  7  x  10“8  seconds  corresponds  to  a  vacuum  pressure  of  2.0  x  10-17  Pa, 
or  1.6  x  10"1®  torr. 
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Ill  Simulation  of  a  Cryogenic  Antihydrogen  Storage  System 

The  analysis  given  in  the  previous  section  provides  the  basis  for  a  computer 
program  to  determine  the  cooling-power  requirements  of  a  simple  antihydrogen 
storage  system.  In  this  program,  several  system  parameters  are  varied  to  determine 
their  effect  on  cooling  power:  the  temperatures  of  the  antimatter  ball  and  the 
chamber  wall,  the  vacuum-chamber  pressure,  the  chamber  radius,  and  the  mass  of 
the  antimatter  ball,  and  the  emissivity  of  solid  hydrogen. 

Some  needed  numbers  are  not  available,  including  the  graybody  emissivity  of 
solid  hydrogen  and  the  precise  spectra  of  the  proton- antiproton  annihilation 
products.  Apparently  solid  hydrogen's  emissivity  is  still  unknown.  In  his  analysis, 
Gaines  assumes  an  emissivity  of  0.5  without  justification.  In  this  program,  the 
emissivity  is  simply  varied  as  a  parameter.  At  this  writing,  precise  figures  have 
not  yet  been  obtained  on  the  energy  spectra  of  proton-antiproton  annihilation 
products.  These  data  exist,  but  do  not  seem  to  be  published  in  numerical  form. 
As  a  rough  estimate,  an  exponential  fit  was  done  to  a  number  of  values  read  from 
the  spectra  in  Gaines'  report.  On  normalizing  this  fit,  it  becomes  apparent  that 
Gaines'  graphs,  which  he  obtained  by  digitizing  an  unsealed  graph  in  Forward's 
report,  give  probabilities  which  are  too  high  by  a  factor  of  ten  (though  Gaines' 
results  based  on  the  spectra  are  correct).  The  fits  used  for  this  study  are 
correctly  normalized. 
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1  Description  of  the  Implementation 

The  program  uses  a  finite-difference  method,  describing  the  system  in  terms 
of  incremental  energy  releases  over  short  time  intervals.  The  energy  incident  on 
the  chamber  wall  is  evaluated  over  a  brief  time  increment  and  divided  by  the  time 
to  determine  the  rate  of  energy  influx  and  thence  the  cooling  power  required 
during  that  time.  The  chosen  time  increment  is  1  ms,  the  time  required  for  a  heat 
impulse  to  subside  to  1%  of  its  peak  value  at  a  distance  of  1.2  cm  from  the  event. 
As  will  be  shown  in  the  next  section,  the  regime  of  operation  in  this  study  was 
effectively  steady-state:  at  the  collision  rates  considered,  the  ice  ball  temperature 
is  constant  over  any  reasonable  system  lifetime.  Therefore,  for  each  set  of 
parameters,  the  results  after  the  first  iteration  are  considered  to  hold  continuously. 
The  program  can  be  divided  in  two  parts:  a  preliminary  portion  where  the  environment 
is  set  up,  and  a  working  portion  in  which  a  single  parameter  is  varied. 

Preliminary  Setup 
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Range  (cm) 


Figure  III- 1 .  Pton  Ranges  in  Iron. 

Before  beginning  to  vary  the  chosen  parameter,  various  environmental  quantities 
must  be  calculated.  The  ice  ball's  radius  is  calculated  from  its  mass,  and  then 
used  in  equations  (5)  and  (11)  above  to  determine  the  annihilation  rate  Raj  and 
the  energy  deposition  per  annihilation  Eai  The  product  of  Raj  and  Eaj  is  the 
annihilation  power  to  the  ice  ball  Pal.  The  charged-pion  energy  deposition  to 
the  chamber  wall  is  computed  by  the  residual-range  method  with  the  pion  range 
function  approximating  published  range  data19  with  the  linear  least-squares  fit 
RANGE  =  (0.082757  cm/MeV)*  -  (3.19505  cm),  as  shown  in  Figure  III— 1.  Gamma  rays 
are  attenuated  exponentially,  using  the  cross  section  curve  shown  in  Figure  II-2: 
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(26) 


\h  -5.4549037  1.9320344  17.519744  1.901761  1 

p"  (Inf)2  +  Inf  +  f2  +  f 

+  (2.851937-  10"5)f-0. 1458238 

where  E  is  in  MeV  and  \i,/p  is  in  cm2/g.  Gaines'  figures,  appropriately  modified 

as  discussed  in  the  previous  chapter,  are  used  to  evaluate  annihilation  energy 
incident  on  the  antimatter  sphere.  In  the  absence  of  numerical  data  on  pion  and 
gamma  energies,  the  spectra  were  approximated  as 

/(f  )=  AeBE  +  CeDE  +  F  (Inf  )2  +  Glnf  +  HE  2  +  IE  +  J  (27) 

where  E  is  in  MeV  and  the  coefficient  values  are  given  in  Table  III-l. 


|  Table  III-l.  Exponential  Pit  Coefficients  for  Pion  and  Gamma  Spectra  | 

Coefficient 

Pion  Spectrum 

Gamma  Spectrum 

A 

-0.1238 

-0.1411 

B 

9.9572  X  10"7 

1.08615  x  10“6 

C 

7.071  x  10~2 

5.578  x  10~2 

D 

-1.0187  x  10"6 

-1.0481  x  10~6 

F 

1.018  x  10~3 

-1.632  x  10-3 

G 

-6.106  x  10"3 

1.443  x  10~2 

H 

2.809  x  10-8 

7.985  x  10~9 

I 

-5.666  X  10"5 

-9.453  x  10-8 

J 

7.867  x  10-2 

5.709  X  10~2 
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Variation  of  Parameters 

For  each  value  of  the  parameter  under  study,  a  single  iteration  of  the 
time -dependent  model  is  run.  This  entails  computing  the  radiant  power  q  using 
equation  (3),  then  adding  it  to  the  annihilation  power  Pa£  to  find  the  total  cooling 
power  Pj  as  shown  inequation  (2).  The  temperature  increase  rate  dTj/dt  is  computed 
from  equation  (14),  with  the  heat  capacity  cvj  given  by  equation  (15).  If  further 
iterations  were  to  be  done,  the  temperature  would  then  be  increased  by  {dTj/dtfdt, 
the  time  incremented  by  dt,  and  the  process  begun  again  with  the  new  temperature 
Tj.  At  the  low  annihilation  rates  considered,  however,  dTj/dt  is  approximately 
zero. 

2  Summary  of  Relevant  Quantities 

The  quantities  dealt  with  can  be  categorized  according  to  whether  they  are 
given  or  assumed,  calculated  (derived),  or  are  variable  parameters.  A  summary  of 
the  relevant  quantities  appears  below. 

Known  Quantities 

The  Debye  temperature  of  solid  hydrogen  120.3/C)20 
Hydrogen  and  iron  cross  sections  for  pions  and  gammas. 

Annihilation-product  energy  spectra. 

Calculated  Quantities 

Specific  heat  of  solid  hydrogen  as  a  function  of  temperature  (Cyi). 
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Antihydrogen  temperature  (Tj)  as  a  function  of  time. 

Radius  of  the  antihydrogen  ball  (mj) 

Annihilation  energy  deposition  rate  in  antihydrogen  (Paj). 

Annihilation  energy  deposition  rate  in  the  wall  (Pa2)- 
Parameters  to  be  Varied 

A  total  of  seven  parameters  are  varied  one  at  a  time.  When  a  parameter  is 
not  being  varied,  it  is  given  a  baseline  value  selected  to  place  the  system  in  a 
near  steady-state  condition.  The  parameters  and  their  baseline  values  are: 

1.  Constant  pressure  of  vacuum  (Pv).  This  parameter  is  varied  from  10~40 
tore  to  1  torr,  with  a  value  of  10-19  torr  in  the  baseline  case. 

2.  Mass  of  the  antihydrogen  ice  ball  (rj).  This  parameter  is  varied  between 
1  ng  and  1  kg,  with  10  mg  as  the  baseline  case. 

3.  Radius  of  the  vacuum  chamber  This  parameter  is  varied  from  1.1  to 

100  times  the  ice  ball  radius,  and  has  a  baseline  value  of  0.1  m. 

4.  Constant  temperature  of  the  chamber  wall  (T2).  This  parameter  is  varied 
from  0.01  K  to  1.50  K.  Its  baseline  value  is  1.00  K. 

5.  Initial  temperature  of  the  ice  ball  (Tqi).  This  temperature  also  runs 
from  0.01  K  to  1.50  K,  but  its  baseline  value  is  1.50  K. 
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6.  Emissivity  of  solid  hydrogen  fe,).  This  quantity  has  apparently  not  been 

■easured,  so  it  is  varied  as  a  parameter  between  0.001  and  1.000,  with  a 
baseline  value  of  0.050. 

7.  Thickness  of  the  chamber  wall.  A  baseline  value  of  1.0  cm  is  used,  but 
this  parameter  is  varied  between  0.5  cm  and  3.0  cm  to  allow  for  different 
vacuum-chamber  designs. 
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IV  Results  of  Simulation 


1  Effect  of  Chamber  Pressure 


1 .0E-071 .0E-061 .0E-051 .0E-041 .0E-031 .0E-021 .0E-011 .0E+001 .0E+01 
Chamber  Pressure  /  le-19  torr 

- Total  Power  - Radiation  Power  .  Annihilation  Power 

Figure  IV-la.  Cooling  Power  Requirement  at  Various  Chamber  Pressures 

The  pressure  to  which  the  vacuum  chamber  can  be  evacuated  Is  the  key 
parameter  of  the  seven  investigated  here.  Figure  IV-la  shows  the  effect  of  chamber 
pressure  (in  units  of  10-19  torr)  on  system  cooling  power.  The  figure  shows  the 
radiative  and  annihilation  components  of  the  cooling  power,  demonstrating  that 
operation  above  about  10-23  torr  is  dominated  by  annihilation,  and  operation  below 
about  10-24  torr  is  dominated  by  thermal  radiation.  The  cooling  power  does  not 
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approach  1  microwatt  until  pressure  Is  above  10'18  torr.  The  following  variations 
of  the  other  parameters  were  run  at  four  chamber  pressures.  In  addition  to  the 
baseline  case  of  10-19  torr,  10-14,  10'16,  and  10-22  torr  were  examined. 


dT/dt  ( 1 e - 40  K/s) 


Chamber  Pressure  (le-19  torr) 

Figure  IV -lb  Temperature  Increase  Rate  at  Various  Chamber  Pressures 

Another  factor  of  interest  is  the  temperature  increase  rate  of  the  ice  ball, 
dTj/dT.  Figure  IV- lb  shows  that  the  system  can  indeed  be  considered  steady-state 
at  the  chamber  pressures  under  consideration  here. 
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2  Effect  of  Antihydrogen  Mass 


Ceding  Power  (W) 


Particle  Mass  (kg) 

■—  le-19  torr  - le-14  torr  . le-16  torr  - le-22  torr 


Figure  IV-2.  Cooling  Power  Requirement  for  Various  Antihydrogen  Masses 


Figure  IV-2  shows  the  effects  of  pressure  and  ice-ball  mass  on  the  system 
cooling  power.  As  one  might  expect,  the  mass  of  antimatter  contained  has  a 
significant  influence  on  system  cooling  power.  There  is  no  qualitative  difference 
in  the  effect  of  mass  at  high  and  low  pressures. 


3  Effect  of  Chamber  Radius 

The  radius  of  the  vacuum  chamber  effects  the  cooling  power  only  when  the 
chamber  radius  is  very  close  to  the  ice  ball  radius  and  the  system  is  operating 
at  a  pressure  such  that  the  radiative  component  of  the  cooling  power  is  significant. 
This  effect  can  be  seen  in  the  10"22-torr  line  in  Figure  IV-3. 


34 


Cooling  Power  (W) 


<“““  le-19  torr  - le-14  torr  le-16  torr  —  le-22  torr 


Figure  IV -3.  Cooling  Power  Requirement  for  Various  chamber  Radii 

4  Effect  of  Wall  Teaperature 

Figure  IV-4  shows  that  with  the  ice  ball  fixed  at  1.5  K,  varying  the  wall 
temperature  produces  little  change  in  cooling  power.  The  effect  is  soaewhat 
greater  if  the  pressure  places  the  sy9tea  in  a  radiation-doainated  power  regiae. 


35 


Cooling  Power  (W) 


— >  le-19  torr  -  le-14  torr  -  le-16  torr  - le-22  torr 


Figure  IV -4.  Cooling  Power  Requirement  at  Various  Wall  Temperatures 


5  Effect  of  Antihydrogen  Temperature 

The  ice  ball's  temperature  is  a  significant  factor  only  at  low  pressures,  when 
radiation  dominates  annihilation  power.  Otherwise,  it  has  no  effect  on  the  cooling 
power,  as  shown  in  Figure  IV-5. 


36 


Cooling  Power  (W) 


Figure  (V-5  Cooling  Power  at  Various  Antihydrogen  Temperatures 


6  Effect  of  Hydrogen  Eaisslvity 

Eaissivity  is  another  factor  that  is  significant  only  at  low  pressures,  as 


shown  in  Figure  IV-6. 


Cooling  Power  (W) 


—  le-19  torr  - le-14  torr  - le-16  torr  - le-22  torr 


Figure  IV -6.  Cooling  Power  Requirement  at  Various  Emissivities 

7  Effect  of  Wall  Thickness 

Wall  thickness  becoaes  sore  important  at  higher  pressures,  when  annihilation 
energy  is  the  dominant  contributor  to  cooling  power.  The  effect  is  not  large  at 
the  pressures  and  thicknesses  considered  here,  however,  as  shown  in  Figure  IV-7. 
It  is  not  apparent  in  the  figure,  but  the  curves  are  linear  at  all  pressures, 
indicating  that  3  cm  is  much  less  than  the  interaction  lengths  of  the  annihilation 
radiation. 
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Cooling  Power  (W) 

I.0E-04 
1 .0L-05 
1 .0E-06 
1 .01-07 
1.0E-08 
1.01-09 
1.0E-10 
J  .01- 1  I 
1.0E-J2 

O.OOS  O.Ol  0.015  0.02  0.025 

Wall  Thickness  (m) 

—  le-19  torr  -  le-14  torr  . le-16  torr  - le-22  torr 

Figure  IV-7.  Cooling  Power  Requirement  for  Various  Wall  Thicknesses 


V  Analysis  and  Conclusion 

It  nay  be  possible  to  store  antimatter  In  the  manner  discussed  In  this  study, 
but  too  many  issues  remain  open  to  make  a  definitive  judgement  one  way  or  the 
other.  Subject  to  the  many  explicit  and  implicit  assumptions  made,  it  may  be 
possible  to  store  useful  quantities  of  antimatter  in  the  manner  discussed  without 
exceeding  the  cooling-power  limits  of  present  cryogenic  techniques.  However,  there 
are  problems  with  the  assumptions,  some  of  which  are  highlighted  below. 

For  the  most  part,  all  of  the  pressure  levels  discussed  are  below  the  measurable 
limit  of  10~14  torr.  It  is  by  no  means  a  given  that  such  pressures  can  be  sustained. 
Even  should  the  required  temperatures  and  number  densities  be  reached,  the  concept 
of  pressure  may  not  be  applicable  in  the  sense  in  which  it  is  used  here:  that  is, 
the  Maxwell-Boltzmann  distribution  and  the  ideal  gas  law  might  not  apply.  The 
question  of  just  where  the  gas  in  the  vacuum  chamber  comes  from  is  also  left  open. 
It  would  be  reasonable  to  think  that  any  hydrogen  gas  in  the  system  got  there  by 
sublimation  from  the  chamber  wall.  Were  this  the  ca„e,  however,  there  would  be  a 
similar  sublimation  rate  from  the  antihydrogen  sphere,  and  disaster  might  be  near. 
Perhaps  also  3He  from  the  cooling  bath  outside  the  chamber  wall  would  be  the  more 
likely  gas  to  appear  in  the  system. 

The  assumption  that  the  antihydrogen  was  entirely  in  the  para  state  has 
significant  Impact  on  the  results,  and  may  turn  out  to  be  ill-founded.  Parahydrogen 
has  very  different  physical  properties  from  orthohydrogen,  including  density,  thermal 
conductivity,  vapor  pressure,  and  phase-change  temperatures  and  enthalpies. 
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Orthohydrogen  also  is  a  significant  source  of  stored  energy  at  low  temperatures, 
since  it  releases  a  large  amount  of  energy  (520  J/g  or  0.011  eV/molecule)  when 
it  converts  to  parahydrogen.  This  is  comparable  to  the  energy  deposited  per 
molecule  by  a  charged  pion  from  a  proton- antiproton  annihilation,  as  shown  by 
equation  (19)  in  Section  II.  If  the  process  by  which  antihydrogen  may  someday  be 
solidified  (a  crucial  assumption  in  itself)  turns  out  to  produce  an  end  product 
with  significant  concentrations  of  anti -orthohydrogen,  this  analysis  would  certainly 
have  to  be  redone  to  include  the  effects  of  the  mixed  substance. 

The  assumptions  that  the  ice  ball  has  a  single  temperature  and  that  no 
sublimations  occur  are  also  suspect.  There  are  several  related  problems  in  this 
area.  First,  more  analysis  needs  to  be  done  on  the  process  of  energy  deposition 
by  charged  pions  in  solid  hydrogen.  The  possibility  remains  that  pions  will  cause 
local  sublimation  of  hydrogen  molecules.  Second,  since  the  Impulse  response  time 
of  the  hydrogen  may  be  comparable  to  the  annihilation  rate,  hot  paths  within  the 
ice  ball  may  overlap.  Since  the  warmer  regions  are  more  conductive,  heat  may  be 
concentrated  at  Impact  points  on  the  surface,  where  sublimation  could  occur  into 
the  chamber.  Third,  the  question  of  shock  effects  on  the  ice  ball  from  annihilations 
has  been  ignored.  Finally,  the  hydrogen  vapor  pressure  curve  used  was  fit  to  data 
points  which  must  themselves  have  been  extrapolations;  its  accuracy  may  be  poor. 

The  cooling-power  requirements  of  storing  solid  antimatter  for  extended  times 
may  not  be  an  obstacle,  if  the  proper  conditions  obtain.  However,  whether  these 


41 


conditions  are  indeed  possible  remains  in  doubt.  Many  difficult  questions  remain 
to  be  answered  before  it  can  be  stated  that  the  goal  of  long-term  storage  is  or 
is  not  within  the  reach  of  present  cryogenic  technology. 
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